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Abstract 

We lay out the foundations of the theory of second-order confor- 
mal superintegrable systems. Such systems are essentially Laplace 
equations on a manifold with an added potential: (A n + V(x))^ = 
0. Distinct families of second-order superintegrable Schrodinger (or 
Helmholtz) systems (A^ + V'(x))^ = can be incorporated into a 
single Laplace equation. There is a deep connection between most of 
the special functions of mathematical physics, these Laplace confor- 
mally superintegrable systems and their conformal symmetry algebras. 
Using the theory of the Laplace systems, we show that the problem of 
classifying all 3D Helmholtz superintegrable systems with nondegener- 
ate potentials, i.e., potentials with a maximal number of independent 
paprameters, can be reduced to the problem of classifying the orbits 
of the nonlinear action of the conformal group on a 10-dimensional 
manifold. 



1 
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1 Introduction 

From our point of view special functions are, in large part, "special" because 
they arise from mathematical models of physical systems that are completely 
solvable analytically and algebraically. Intuitively we consider such systems 
to be of high symmetry, but this symmetry may be "hidden", i.e., not obvi- 
ous. Special function theory can be based on the notion of superintegrability; 
it is the best concept to date to capture both hidden symmetry and just those 
systems whose associated functions are interesting and useful enough to be 
considered "special" . An n-dimensional Hamiltonian system (2n-dimensional 
phase space), classical or quantum, is integrable if it admits n functionally 
independent commuting constants of the motion (or symmetry operators), 
polynomial in the momenta. It is superintegrable if it is integrable and ad- 
mits 2n-l constants of the motion (the maximum possible but, of course, 
not all commuting). If the functionally independent constants of the motion 
can all be chosen of order k or less in the momenta (or in the derivatives 
for the quantum case) the system is called kth order superintegrable. Su- 
perintegrability is a much stronger requirement than integrability, indeed 
superintegrable systems can be solved algebraically. 

Special functions are connected to superintegrable systems in several 
ways. The most obvious is that they occur when one computes the eigenfunc- 
tions of the Hamiltonian operator or other symmetry operator in a 1st or 2nd 
order quantum superintegrable system. The superintegrability forces variable 
separability, usually in multiple systems. The majority of special functions 
of mathematical physics arise from separable coordinates in this way. The 
symmetry algebras generated by 1st order quantum superintegrable systems 
are Lie algebras and the typical special functions arising are spherical har- 
monics, and other orthogonal polynomials (22j [2Tj . Superintegrable systems 
of 2nd order are multi-integrable (so multi-separable) and distinct classes of 
special functions can be related to one another in a single system. Indeed, 
the basic properties of Gaussian hypergeometric functions and their various 
limit ine; CclSGS, clS well as Lame, Mathieu and Heun functions, and ellipsoidal 
harmonics are associated with 2nd order superintegrable quantum systems 
via separation of variables. The algebra formed by the generating symmetries 
again closes under commutation. For example, consider operator superinte- 
grable systems of the form (A + V)^f = E^> on a 2-dimensional conformally 
flat manifold with potential function V, and symmetry generators of order 
no more than 2. Then the system of symmetries closes, sometimes at order 
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1 (Lie algebras and Lie groups, V = 0), sometimes at order 3 for degenerate 
potentials, and sometimes at order 6 for non-degenerate potentials. There 
are no other possibilities. Each such system is multiseparable. (For n > 2 
similar statements hold but there are more possibilities). Closure at order 
1 corresponds to a Lie algebra. The monograph [19], written before the 
word "superintegrable" was coined, is really about some simple 2nd order 
super integrable systems whose algebra closes at order 1. 

Closure at orders 3 or 6 defines quadratic algebras (NOT Lie algebras) 
whose algebraic representation theory gives crucial information about the 
possible energy eigenvalues E and the expansion of one integrable eigenba- 
sis in terms of another (i.e., the expansion of one class of special functions 
in terms of another). The representation theory of these quadratic algebras 
is of great intrinsic interest and leads to another connection with the the- 
ory of special functions. For example, consider the system (A + a/(si) 2 + 
b/{s2) 2 + c/(s3) 2 )\l/ = where A is the Laplace-Beltrami operator on the 
2-sphere (si) 2 + (s 2 ) 2 + (S3) 2 = 1. This is 2nd order superintegrable with a 
quadratic algebra that closes at order 6 (i.e., differential operators of order 
6). It is an amazing fact that for eigenvalues E of finite multiplicity, this 
algebra is precisely the structure algebra for the Racah polynomials in their 
full generality. The algebra for the infinite dimensional bounded below rep- 
resentations of the quadratic algebra yields the Wilson polynomials in their 
full generality [TBI E El ED]- Thus all of the classical discrete orthogonal 
polynomials and their Wilson polynomial generalizations appear naturally 
in the representation theory of the quadratic algebra. The interpretation 
here is that the Schrodinger equation separates in spherical coordinates and 
elliptic coordinates. The spherical coordinate eigenf unctions are products of 
Jacobi polynomials, known as Karlin-McGregor 2-variable orthogonal poly- 
nomials in this case. The elliptic coordinate eigenf unctions are products of 
Heun polynomials. The expansion of the Heun polynomial basis in terms 
of the Karlin-McGregor basis can be computed directly from the represen- 
tation theory of the quadratic algebra, but is complicated. One can also 
choose spherical coordinates based on rotation about the 3-axis or spherical 
coordinates based on rotation about the 2-axis. The coefficients giving the 
expansion of one spherical coordinate basis in terms of the other are just the 
Racah (or Wilson) polynomials. Thus, study of the representation theory of 
the quadratic algebra leads to classes of special functions, such as orthogonal 
polynomials of a discrete variable, that do not arise directly from variable 
separation. 

Special functions are also associated with higher order superintegable 
systems. Thus the Painleve transcendents (not associated with variable sep- 
arability) appear in the study of 3rd order superintegrable systems and the 
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representations of their cubic symmetry algebras, El HE]- Similarly, su- 
perintegrability can be related to generalized hypergeometric functions and 
many variable hypergeometric functions (Appell functions, Lauricella func- 
tions, Horn functions, etc.) For 1st order superintegrability some of this is 
in [19], but this is mostly virgin territory for superintegrability of order > 2. 

Superintegrabilty can also be studied for equations of the form (A n + 
V)^ = on conformally flat n-dimensional manifolds, and that is the sub- 
ject of this paper. For example, Volkmer [23] treated the equation (A n + 
a>i/(xi) 2 + ••■ + On/ (x n ) 2 )^ = on n-dimensional flat space, in a study of 
generalizations of ellipsoidal wave functions. Such functions arise when this 
equation is separated in conical coordinates. From our point of view, this is 
a conformally 2nd order superintegrable system with degenerate potential. 
The added insight provided by superintegrability theory is that the conformal 
symmetry algebra of this equation closes at order 6. Thus the representation 
theory of this quadratic algebra will provide important spectral information 
about the symmetry operators of the Laplace equation. 

This paper inaugurates the study of second-order conformally superinte- 
grable operator systems on conformally flat manifolds, essentially Laplace 
equations on a manifold with an added (very special) potential: (A n + 
V(x))^/ = 0. There are several important features of this approach. First, 
distinct families of second-order superintegrable Schrodinger (or Helmholtz) 
systems (A' n + V (x.))^ = Ety can be incorporated into a single Laplace equa- 
tion, and we can exploit the relationship between them. Second, via a gauge 
transformation we can always transform the Laplace problem to flat space 
and make direct use of the conformal symmetry algebra so(n + 2, C) of the 
Laplacian. Using this approach, we will show that the problem of classify- 
ing all 3D Helmholtz superintegrable systems with nondegenerate potentials, 
i.e., potentials with a maximal number of independent paprameters, can be 
reduced to the problem of classifying the orbits of the nonlinear action of the 
conformal group on a 10-dimensional manifold. Eventually, this should lead 
to a new classification structure for special functions and their properties. 

2 The classical Laplace-type conformal super- 
integrable system 

We start by defining Laplace-type conformal superintegrability in classical 
mechanics. The Hamiltonian system is 7i = where 7i = 7{q + V^(x) and 
Ho = Y^j=i 5 ,u ( x )PiPi is the free particle Hamiltonian on a real or complex 
conformally flat pseudo-Riemannian space. The phase space is 2n dimen- 
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sional with local coordinates 



(x,p) = (x ir -- ,x n ,p 1 ,--- ,p n ). 



The Poisson bracket of functions /, g on the phase space is 



n 



if, 9} = ^2(d x Jd Pi g - d Pi fd Xi g) 



i=l 



The condition Ti = restricts us to a 2n — 1 dimensional hypersurface in 
phase space. A conformal symmetry of this equation is a function «S(x, p) 
such that {S,Ti} = Rs(x.,p)Ti for some function R$. Two conformal sym- 
metries S, S' are identified if S = S' + RTi for R any function on phase 
space, since they agree on the hypersurface Ti = 0. The system is confor- 
mant/ superintegrable if there are 2n — 1 functionally independent conformal 
symmetries, <Si, • • • , iS>2n-i with Si — Ti which firstly, are polynomial in the 
momenta and secondly, the symmetries £2, • • • , S^n-i are still functionally 
independent on restriction to the hypersurface Ti — 0. The system is second 
order conformally superintegrable if each of the basis symmetries Si can be 
chosen as a second order polynomial in the momenta. The condition that S 
is a conformal symmetry implies dS/dt = |{«S, H} = \RsH so dS/dt = 
at any point (x, p) on the hypersurface Ti = and S is constant along any 
trajectory satisfying Hamilton's equations 



Note that if a point of the trajectory lies on the hypersurface Ti = then all 
points on the trajectory lie on this hypersurface. Thus for constants c = (q), 
with C] = we can solve the equations <Sj(x, p) = q, i — 1, • • • , 2n — 1 
analytically to get a 1-parameter trajectory. 



3 The operator Laplace-type conformal su- 
perintegrable system 



Here A n is the Laplace-Beltrami operator on a real or complex conformally 
flat Riemannian or pseudo Riemannian manifold. A conformal symmetry 
of this equation is a partial differential operator S in the variables x = 



2«9 t x = d p H, 2d t p 



-dji. 



Systems of Laplace type are of the form 



HV = A n * + = 0. 
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(xi, ■ ■ ■ ,x n ) such that [S, H ] = SH — HS = RsH for some differential 
operator R$. A conformal symmetry maps any solution \1/ of ([T]) to another 
solution. Two conformal symmetries S, S' are identified if S = S' + RH for 
some differential operator R, since they agree on the solution space of (0Q). 
The system is conformally superintegrable if there are 2n — 1 functionally 
independent conformal symmetries, Si, - • - , S^n-i with Si = H. It is second 
order conformally superintegrable if each symmetry Si can be chosen to be 
a differential operator of at most second order. 

We can distinguish three types of conformally superintegrable Laplace 
equations. The first type is just a recasting of a Helmholtz superintegrable 
system H'^> = E^> into Laplace form H^f = where H = H' — E. Thus the 
parameter E is absorbed into the potential. The second type is a restriction 
of a Helmholtz superintegrable system H'ty = Ety where H' = A n + V to a 
fixed energy eigenvalue E Q The resulting system 

(A n + V(x) -E )V = 

is trivially conformally superintegrable. Here E is a fixed constant in the 
potential V = V — E , not a parameter, and it is not permitted to add 
further nonzero constants to the potential, [1] . Generically, all of the sym- 
metries of a system of this type will be those inherited from the Helmholtz 
equation, so the restriction seems to be no special interest. However, for 
some particular energies Eq new truly conformal symmetries may appear so 
that the structure of the symmetry algebra will change. The third type of 
conformal superintegrable system is one of the form ([T]) where the Helmholtz 
equation if \I/ = E^ is not superintegrable. In this case the truly conformal 
symmetries necessarily appear. 

3.1 A Laplace equation example with degenerate po- 
tential 

For our first example we consider the equation 

H * - XX + §)* = °- ( 2 ) 

i=l * 

(This equation was treated by Volkmer in his study of generalized ellipsoidal 
harmonics in n dimensional Euclidean space, [23].) In this case the Helmholtz 
equation H^> = E^> is second order superintegrable with degenerate poten- 
tial. A basis of generators for the quadratic algebra of symmetries is given 
by the n(n + l)/2 second order symmetries 

P i = % i + %, J = l,---,n (3) 
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and 

x 2 x 2 

J jk = (xjd Xh - x k d X] f + a,— 1 + a&-f , 1 < j < k < n. (4) 

X j X k 

(Of course there are functional relations between these symmetries since only 
a 2n — 1 element subset is functionally independent. For the case n = 3 these 
relations can be found in [T7], by restriction.) What makes this potential of 
particular interest is that for E = the system admits new (truly conformal) 
symmetries. The most obvious is the dilation symmetry 

D = -J>^,-^-. (5) 

(The constant term has been added for convenience in the calculations to 
follow.) Note that [D, H] = 2H. Further, there is a non-local symmetry / 
defined by 

1 n 
"<*,», = r" = J>?. (6) 

1=1 

Here, 

[J,tf]=r 4 tf, / = /-!, 

so if 5 is a differential symmetry, so is ISI~ l . Now IJj^I' 1 = Jj k and 
IDI^ 1 = —D, so we get nothing new. However the operators Kj = IPjI~ l 
are new conformal symmetries: 



— n ^2\2 



^ = (n - 2) Xj - x % 3 + 2x i Y Xid ^ + a i i= \ Xi 



X 

i=l i=l / 3 



(7) 



These symmetries are not independent of one another. In particular we have 
the identities 

n n n 2 

5^ = 0,5^ = 0, y Jik+D 2 + Y^-^- = n, (s) 

i=l i=l l<i<fc<" *=1 

each valid on the solution space of H^ff = 0. The first order conformal 
symmetry D acts on the second order symmetries via 

[D, P,\ = 2P V [D, K s \ = -2K V [D, J jk ) = 0. (9) 

We also have the second order commutator relations 

[P u Pj\ = 0, [K u K J ] = 0. 



7 



The expressions for the commutators [Pi, Jjk], [Ki, Jjk] and [Pi, Kj] are more 
complicated and the structure of the symmetry algebra generated via com- 
mutation is not completely clear at this time. Note however, since (CQ) 
can be thought of as a restriction of the singular isotropic oscillator for 
the Helmholtz equation in n dimensions, the operators [Pi, Jjk][Pi>, Jj'k'] + 
[Pi>, Jfk'][Pi, Jjk] can De expressed as symmetrized third order polynomials 
in Pi and Jhm- Similarly by applying the symmetry / we see that the opera- 
tors [Ki, Jjk)[Ki>, Jj'k'] + [Kir, Jj'k'][Ki, Jjk] can be expressed as symmetrized 
third order polynomials in Ki and Jhm- Using the same reasoning we see 
that fourth order operators of the form [[Pi, Jjk], Jj'k'] or [[Pi, Jjk], Pi'] can be 
expressed as symmetrized second order polynomials in Jj» t k" and Pj«. Simi- 
larly fourth order operators of the form [[K, Jj k ], Jj'k'] or [[Ki, Jjk], K^] can 
be expressed as symmetrized second order polynomials in Jj"^k" and K^r. 

It is important to recognize that for eigenfunctions ty\ of the dilation 
operator, D^/\ = \^\, the conformally superintegrable system on flat space 
specializes to the Helmholtz equation on the n — 1 sphere with generic po- 
tential: 

J^x = (-\ 2 + ^)*x, (10) 

l<j<fe<71 

a superintegrable system. Thus the conformal symmetry algebra of (T5]) can 
be regarded as a dynamical symmetry algebra for ( flOl) since in general these 
conformal symmetries will change the eigenvalue A, hence the energy. Since 
the operators [Pi,Kj] commute with D, they are symmetries of (imp hence 
expressible in terms of the commutators of the basis symmetries Jjk for the 
n — 1 sphere with generic potential. Similarly, since symmetries of the form 
[Pi, Jjk][Ki', Jj'k'] + [K', Jj'k'][Pi, Jjk] also commute with D they must be ex- 
pressible as third order symmetrized polynomials in the basis symmetries 
Jjk- In like manner, fourth order symmetries for the form [[Pi, Jjk],K(\ and 
any like expressions that commute with D must also be expressible as sym- 
metrized second order polynomials in the Jj'k 1 - 

Putting these observations together, we conclude that the conformal sym- 
metry algebra of (j2J) generated by Pi, Jjk, Kp, D must close at order six, so 
it is a true quadratic conformal symmetry algebra. The structure theory for 
the quadratic conformal symmetry algebra and the complete set of functional 
relations among the generators is yet to be determined. However, for n = 2, 3 
we can understand the functional relations. 

For n — 2, the simplest and atypical case, there are 3 functionally inde- 
pendent generators, whereas we have 6 second order conformal symmetries 
P U P 2 , Ji 2 , K U K 2 ,D 2 . The relations 

Pi + P 2 = H ~ 0, K 1 + K 2 = (xl + x\fE ~ 0, 
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Ji 2 + D 2 + ai + a 2 = (x 2 + a^)# ~ 0, 
^12 " ^(KiPi + Pi^i) - 5Ji 2 - 3(ai + a 2 ) - 4 ai a 2 = (11) 
- (x 2 (x 2 - 2y 2 )dl - x'd 2 y + 4x 3 yd xy 

+x(6a; 2 - 4?/ 2 )9 x + 6x 2 yd v + 9a; 2 + a x x 2 y 2 + 2aiy 2 - a 2 — ) H ~ 

yield the complete structure, where we write A ~ P if the operators A, 5 have 
the same action on the null space of H. Indeed, we can take H, Pi, K 2 as the 
functionally independent generators. Then we have commutation relations 
(ED and 

[Pi, K x ] ~ P 3 + 4(1 + 2a! + 2a 2 )P 

which determine everything. In this case the algebra closes at order 3. 

If we further set a 2 = then the structure of the conformal symmetry 
algebra changes again. We get a new first order symmetry L\ = d X2 and, since 
K 2 = L\ becomes a perfect square, another first order conformal symmetry 

L2 = ~{x\ + x 2 2 )d X2 + 2x 2 (x 1 d Xl + x 2 d X2 ). 

Thus the system is now first order superintegrable with structure 

[D,Lx] = Lu [D,L 2 ] = -L 2 , [L U L 2 ] = -2D, 

the Lie algebra s£(2). To within a gauge transformation the Laplace equation 
= is just a complexification of the EPD equation, studied in [13] from 
the group theoretic point of view. 

For n = 3 there are 10 second order symmetries Pi, P 2 , P3, K\, K 2 , K%, 
J12, J13, J23, D 2 , only 5 of which are functionally independent. This is ex- 
plained by the 3 second order relations (jSJ) and 2 eighth order relations, one 
relating the Pj, Jjk and one relating the K iy J^. The algebra closes at order 
6, as it does for all n > 2. 

3.2 Laplace equation examples with nondegenerate 
potential 

Example 1 We choose three dimensional flat space with Cartesian coordi- 
nates x, y, z, though the construction works in all higher dimensions. We 
set 

— Sl — S2 — 53 2 1 2 , 2 1 2 _ -1 

1 + s 4 1 + s 4 1 + s 4 
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essentially the stereographic projection of the three-sphere on Euclidian space. 
We now choose generic coordinates on the three-sphere that correspond to R- 
separable coordinates for the flat space Laplace equation, as in J7P) /: 



{H - a)(y - a){p - a) Up, - b){u - b)(p - b) 

S| W (b-a){a-l)a ' * 2 V (a -&)(&- 1)6 : 




(p - l)(u - \){p - 1) jpiTp 

S4 -y-^b- (12) 



= dx 2 + dy 2 + dz 2 = - - ^ 2 (ds\ + ds\ + ds\ + ds 2 ) . (13) 
We can invert the coordinate transformation according to 



2x 2y 2z 

y 2 + z 2 + V 

1 — x 2 — y 2 — z 2 



1 x 2 + y 2 + z 2 + 1 ' 2 x 2 + y 2 + z 2 + 1 ' 3 x 2 + y 2 + z 2 + 1 



4 x 2 + y 2 + z 2 + 1 

If we proceed at the classical level, then by choosing the coordinates indi- 
cated above we can achieve multiseparation and hence conformal superinte- 
grability with a potential V of the form 

V = (l + stf(% + % + % + %-a B ) (14) 
\s 1 s 2 s 3 s 4 / 

0,1 _|_ °2l _|_ Hi _|_ 4 ° 4 4fl5 



x 2 y 2 z 2 (1 — x 2 — y 2 — z 2 ) 2 (1 + x 2 + y 2 + z 2 ) 2 

This is a conformally superintegrable system with nondegenerate potential. 
For the quantum analogue we obtain a quite specific spectrum of the corre- 
sponding Laplace operator on the three sphere with eigenvalue shifted by 3/4. 
Indeed, writing Laplace's equation 

(d xx + d yy + d zz + V)^ = (15) 

in terms of the ellipsoidal coordinates fi, u, p and applying the gauge transfor- 
mation = (1 + S4)~ 1 / 2 $ we find that ( 1731) transforms to a separable system 
with separation equations 

i . ... .... -i \ a ( o2 1/ 1 1 1 1,„ (a — b)(a — l)a 

(A- a )(A-&)(A-l)A ^ + + _ + _ + - a / ^ 
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(b — a)(b— l)b (1 — b)(l — a) a^ab , 3., 9 , \ 
-a 2 - A _ fe ; - a 3 ± ^ + - (a 5 + -)A 2 + «iA + k 2 J A(A) = 

/or X = p,u, p and separation constants Kx,k 2 - In terms of the function $ 
i/ie Laplace equation becomes 



a* + £ + 3 + 3 + 3-(«>. + §»» = o (16) 

*1 b 2 6 3 6 4 ^ 



where As 3 is the Laplace- Beltrami operator on the three sphere. This clearly 
illustrates the connection with the superintegrable system with nondegenerate 
potential on the three sphere for the energy translated by the specific value 3/4. 
Here the 3/4 is due to the quantization procedure. In 2 dimensions it it doesn't 
appear and for n > 2 dimensions it is 3/16 of the scalar curvature of the 
manifold. The quadratic algebra structure for the system [TB) is preserved by 
the gauge transformation to give the conformal quadratic algebra for system 
173)) . This map from a nondegenerate Helmholtz superintegrable system on 
the 3-sphere with potential 

V = — + — + — + — -a 5 

S l S 2 S 3 S 4 

to a nondegenerate Laplace conformally superintegrable system in flat space 
with potential 

V = {l + s 4 fV = {l + s 4 f (^ + ^ + ^ + ^-a 5 

\s 1 s 2 s 3 s 4 

is not a Stdckel transform, |3J \3i- (It is an example of a Stackel multiplier, 
fT5f.) However, since (1 + S4) 2 is a specialization of the Laplace potential, 
the inverse transform can be considered a Stackel transform. Further, we 
can apply Stackel transforms to this example, and to our previous example, 
to obtain multiple conformal superintegrable systems on manifolds with non- 
constant curvature. This idea will be used to characterize arbitrary Helmholtz 
superintegrable systems on conformally flat spaces as Stackel transforms of 
flat space Laplace superintegrable systems. 

Example 2 A second approach to example (fl^P uses general cyclidic coor- 
dinates p, p, v . We choose 

„2 _ (p-e 1 )(/i-e 1 )(z/-ei) 2 _ {p - e 2 ) (p - e 2 ) (u - e 2 ) 

x 2 



(ei - e 2 ) (ej - e 3 ) (ex - e 4 ) (e x - e 5 ) ' (e 2 - e l )(e 2 - e 3 ) (e 2 - e 4 ) (e 2 - e 5 ) 

(P - e 3) (ji ~ e 3 ) (1/ - e 3 ) 2 (p - e 4 ) (/1 - e 4 ) (1/ - e 4 ) 



(e 3 - e 2 )(e 3 - ei)(e 3 - e 4 )(e 3 - e 5 )' X4 (e 4 - e 2 )(e 4 - e 3 )(e 4 - ei)(e 4 - e 5 ) 
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(p - e 5 )(/i - e 5 )(u - e 5 ) 



Here e±, ■ 
the cone 



' h (e 5 - e 2 )(e 5 - e 3 )(e 5 - e 4 )(e 5 - ei)' 
es are constants. These are the pentaspherical coordinates on 

(17) 



2 i 2 i 2 i 2 i 2 /~i 



and they can be written in terms of projective coordinates X, Y, Z, T 
xi = 2XT, x 2 = 2YT, x 3 = 2ZT, x 4 = X 2 + Y 2 + Z 2 - T 2 , 

x 5 = t(X 2 + Y 2 + Z 2 + T 2 ). 



The Cartesian coordinates x, y, z are given by 



X 

x = — 
T 



X 4 + tx 5 
We also note the relations 



x 2 + y 2 + z 2 - 1 



Y 

y rj-\ 



x 2 



x 3 



x 4 + ix 5 ' 



x 4 + ix 5 



2x 4 



X4 + ix§ 

The metric distance in Euclidean space is 



ds 2 = dx 2 +dy 2 +dz 2 = (x 4 +ix 5 ) 2 



x 2 + y 2 + z 2 + 1 



2ix* 



£4 + ix$ 



(19) 



(20) 



(P - P){p ~ ")dp 2 



(p - u)(n - p)dp 2 



(p - ei)(p - e 2 )(p - e 3 )(p - e 4 )(p - e 5 ) 

(21) 

(z/ — p)(y — p)dv 2 



+ 



(ji - ei)(/i - e 2 )(/i - e 3 )(p - e 4 )(p - e 5 ) (u - ei)(v - e 2 )(^ - e 3 )(v - e 4 )(u - e 5 ) 
VKe can now construct a general potential which is conformally superintegrable 



viz 



V^= (X 4 + 2X 5 ) (-2+-2+-2+-2+Z2) 



2 '2. 
4 x 5 



'1 ^2 ^3 

which is identical to (TJ) when written in terms of Cartesian coordinates. 

Writing Laplace's equation ( T73j) m terms of these coordinates and setting 
\P = (x 4 + ix5) _1//2 $ we obtain a partial differential equation for $ i/iai 
separable: $(p, p, z/) = Ai(p)A 2 (p)A 3 (z/). T/iis /eads to i/ie i/iree separation 
equations 

[4(A - e x )(A - e 2 )(A - e 3 )(A - e 4 )(A - e B )[%+ 

1 ■ 1 -)«J 



k 1 



1 1 

+ 



2 A — ei A — e 2 A — e 3 A — e 4 A — es ' 
(e 4 - e 2 )(ei - e 3 )(e 4 - e 4 )(ei - e 5 ) (e 2 - e 4 )(e 2 - e 3 )(e 2 - e 4 )(e 2 - e 5 ) 

+0l 7T : hCt 2 



(A -ex) 



(A - e 2 ) 
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(e 3 - e 2 )(e 3 - e 1 )(e 3 - e 4 )(e 3 - e 5 ) (e 4 - e 2 )(e 4 - e 3 )(e 4 - ei)(e 4 - e 5 ) 

+o 3 7T v Ha 4 — r 

(A - e 3 ) (A - e 4 ) 

(e 5 - e 2 )(e 5 - e 3 )(e 5 - e 4 )(e 5 - e 4 ) 5 3 

+«5 77 x 7 A + 

(A - e 5 ) 4 



3 

-(ei + e 2 + e 3 + e 4 + e 5 )A 2 + KiA + k 2 



A(A) = 0, 



A = p,/i, 

// we take e^, —> oc then we recover the system above associated with 
the three dimensional sphere. In particular we note that the expression for 
the gauge factor is readily computed. This is a consequence of the relation 
t\x\ + e 2 x\ + e 3 Xg + e±x\ + e 5 xf = — 1, from which we deduce that 

/ , • \2 2 , 2 i 2 , /i 2 2 2\2 /i , 2 , 2 , 2\2 

-(x 4 + ^ 5 ) = e x x +e 2 y +e 3 z + e 4 (l-x -?/ -z ) -e 5 (l + x +y +z ) . 

Example 3 We consider the nondegenerate Laplace conformally superinte- 
grable system in flat space with potential 

ai a 2 z a 3 (x 2 + y 2 - 3z 2 ) a 4 
V(x, y, z) = - — — — + - — — — + — — — + 



(x + iy) 2 (x + iy) 3 (x + iy) 4 (1 — x 2 — y 2 — z 2 ) 2 

(22) 

+ (l+x 2 + y 2 + z 2 ) 2 ' 

(This system will be real in real Minkowski space with coordinates X = x, 
Y — iy, Z = z.) A specialization of this potential is 1 / '(1 + x 2 + y 2 + z 2 ) 2 and 
if we use it to perform a Stackel transform we get a Helmholtz superintegrable 
system on the complex 3-sphere, in complete analogy with our first example. 
This is the system IV in [10]. Another specialization of potential fifty) is 
l/(x + iy) 2 . If we use it to perform a Stackel transform we get a Helmholtz 
superintegrable system on complex flat space, the system IV in FTUtf . 

Example fifty) is just one of a class of nondegenerate supperintegrable 
systems on the 3-sphere and flat space that are induced from nondegenerate 
superintegrable systems on the 2-sphere. Indeed, any superintegrable system 
on the 2-sphere can be imbedded in 3 dimensional flat space in an obvious 
manner. For each 3-parameter potential of a 2-sphere nondegenerate super- 
integrable system, such as listed in [13^ . we get a flat space superintegrable 
system whose coefficients 0/01,02,03 can be read off from the list in fT3^ . 
Then we add the terms a 4 / (1 — x 2 — y 2 — z 2 ) 2 + 05/(1 + x 2 + y 2 + z 2 ) 2 to get 
the corresponding Laplace nondegenerate superintegrable system. 
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3.3 Pentaspherical coordinates and linearization of the 
conformal group 

Here we examine the role of pentaspherical coordinates ([T7]) . (TTSj) . (TH1 . (T2U|) 
in more detail. From these relations we can write 

d x = 2Td Xl + 2Xd X4 + 2iXd xs , d Y = 2Td Xl + 2Yd X4 + 2iYd X5 , 

d z = 2Td xl + 2Zd X4 + 2iZd XB . 

We also recognize d x = Tdx , d y = Tdy , d z = Tdz- From these observa- 
tions we deduce that the spatial derivatives are related to the pentaspherical 
derivatives via 

d x = -(au + ix 5 )d xl + x x (d X4 + id X5 ), d y = -(x 4 + ix 5 )d X2 + x 2 {d X4 + id X5 ), 

d z = -(x 4 + ix 5 )d X3 + x 3 (d X4 + id X5 ). 
The classical analogs of these relations are 

Px = ~(x 4 + ix 5 )p Xl + x x {p XA + ip X5 ), p y = -(x 4 + ix 5 )p X2 + x 2 {p XA + ip X5 ), 

Pz = ~{xi + ix h )p X3 + x 3 (p X4 + ip X5 ). 
From these relations we determine that the flat space free Hamiltonian is 

Ho = P 2 X + p\ + P 2 Z = (xa + ix 5 ) 2 (p 2 xi + p 2 X2 + p 2 X3 + p 2 Xi + p 2 XB ). (23) 

Recalling that motion is restricted to the null cone x\ + x 2 + x\ + x\ + x\ = 0, 
we can consider 7i as a Hamiltonian in 10- dimensional pentaspherical phase 
space with Poisson bracket {-, -}p with respect to which Xi,p Xi are canonically 
conjugate variables only if 

5 5 

{^2 X l> ^0}P = 2(X 4 + X 5 f ^ X kPx k = 0, 
fc=l k=l 

so we restrict our consideration to the subspace of pentaspherical phase space 
that is on the null cone and for which 

xip xi + x 2 Px 2 + x 3 p X3 + x 4 p X4 + x 5 p X5 = 0. (24) 

The point of all this is that the action of conformal symmetries is linearized 
in pentaspherical phase space. Indeed, at the conformal Killing vector level 
we have 

Px = (xip X4 -x A p Xx ) + i(x l p X6 -x 5 p Xx ), p y = (x 2 p X4 -x^p X2 ) + i(x 2 p X5 -x 5 p X2 ), 
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Pz = (x 3 p X4 - x A p X3 ) + i(x 3 p X5 - x 5 p X3 ), 

X Py ~ VPx — %lPx 2 ~ x 2Px\ , UPz — Zp y = X2P XZ — %3Px 2 j 

zp x - xp z = x 3 p Xl - X!p X3 , xp x + yp y + zp z = i(x A p X3 - x 5 p X4 ), 

-(x 2 + y 2 + z 2 )p x + 2x(xp x + yp y + zp z ) = (x x p XA - x±p Xi ) -i(x x p X6 -x 5 p Xl ), 

-{x 2 +y 2 + z 2 )p y + 2y(xp x + yp y + zp z ) = (x 2 p X4 - x A p X2 ) - i(x 2 p Xs - x 5 p X2 ), 

-(x 2 + y 2 + z 2 )p z + 2z(xp x + yp y + zp z ) = (x 3 p X4 - x 4 p X3 ) - i(x 3 p X5 - x 5 p X3 ). 

This means that the classical conformal Killing vectors of 7{q are just those 
of the complex Lie algebra of 5*0(5,0), and the corresponding connected 
component to the identity of the group symmetries is just this linear Lie 
group. Classically, the inversion operation in a sphere corresponds to the 
reflection 

I : xj -> xj, p Xj -> p Xj , j ^ 4, x A -»• -x 4 . p X4 -> -p XA (25) 

Thus the full conformal group action in pentaspherical space is just the linear 
action of 0(5, C). 

A straightforward computation shows that under the null cone restriction 
(JT7J) and the restriction ( J24l) the canonical one-form uo = p x dx + p y dy + p z dz 
on our original phase space goes to the canonical one-form on the restricted 
pentaspherical space: 

5 

uj = p x dx + Pydy + p z dz = '^^p Xk dxk- 

k=l 

Thus 

5 

du = dp x A dx + dp y A dy + dp z Adz = dp Xk A dx^ 

k=l 

so the symplectic two-forms agree and we have achieved an embedding of 
our 6-dimensional phase space into the 10- dimensional pentaspherical phase 
space that preserves Poisson bracket relations. 

The flat space classical system 7i = pi + p 2 + p 2 + V(x) = can now be 
lifted to pentaspherical space: 

5 

H = (x A + zx 5 ) 2 C£pl k ) + V(x)=0 
k=i 
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where now V(x) = (x4 + ix 5 ) 2 V(x) is expressed in pentaspherical coordinates 
Xk- Thus we have the equation 



fc=i 

and the possibilities for V to correspond to a superintegrable system relate 
to properties of confocal quadratic forms in 5-space. This construction was 
exploited long ago by Bocher in his monograph on /^-separation of variables 
for Laplace equations [2]. The general cyclidic coordinates and the coordi- 
nates on the sphere given previously are easily related to this construction. 
To see this observe that general cyclidic coordinates in 5-space are given by 

2 _ (p-ei)(i^-e 1 )(/i-e 1 ) 

X-i 



(ei - e 2 )(ei - e 3 )(ei - e 4 )(ei - e 5 ) 

with analogous expressions for x 2 , ■ ■ ■ ,x$. Now take e\ — and substitute 
A — > 1/A, A = p, p, ^ and ej — > z = 2, 3, 4, 5. We obtain 

2 P V P 2 2 (p ~ e 2 )(v - e 2 )(p - e 2 ) 2 2 (p - e 3 ) (i/ - e 3 ) (p - e 3 ) 
1 e 2 e 3 e4e5' 2 1 (e 2 - e 3 )(e 2 - e 4 )(e 2 - e 5 ) ' 3 1 (e 3 - e 2 )(e 3 - e 4 )(e 3 - e 5 ) ' 

(p - e 4 )(z/ - e 4 )(p - e 4 ) 2 2 (p - e 5 )(i/ - e 5 )(p - e 5 ) 



! ^(e 4 - e 2 )(e 4 - e 3 )(e 4 - e 5 )' 5 1 (e 5 - e 3 )(e 5 - e 4 )(e 5 - e 2 ) ' 

From this calculation we could just as well choose new pentaspherical coor- 
dinates 

X\ = — 1, -X" 2 = ~~ 77 TT ^T' X- 



Xa — 7~~~ TT TT~ ~"T i X, 



2 _ (p-e 2 )(^-e 2 )(p-e 2 ) ^ 2 = (p - e 3 )(z/ - e 3 )(p - e 3 ) 
e 2 - e 3 ) (e 2 - e 4 ) (e 2 - e 5 ) ' 3 (e 3 - e 2 ) (e 3 - e 4 ) (e 3 - e 5 ) 

(p - e 4 ) (i/ - e 4 ) (p - e 4 ) 2 (p - e 5 ) (i/ - e 5 ) (p - e 5 ) 



" 4 (e 4 - e 3 )(e 4 - e 2 )(e 4 - e 5 )' 5 (e 5 - e 3 )(e 5 - e 4 )(e 5 - e 2 ) ' 

from which it is clear that X\ + X 2 + X\ + X| = 1 . If we now relabel these 
coordinates according to X 4 = y 3 , X 5 = y 4 , X 3 = y 2 , X 2 — y 1: X x = y 5 and 
choose Cartesian-like coordinates according to 

V\ 2/2 2/3 
X = ; , y = ; , Z = ; , 

2/4 + W5 2/4 + «2/5 2/4 + W5 

we obtain the separable coordinate systems on the sphere displayed previ- 
ously. What this demonstrates is that the Cartesian-like coordinates are 
determined only by the ratios of the pentaspherical coordinates as indicated 
above. In addition the conformal symmetry group is available, acting on the 
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vector x = (x±, x 2 , x 3 , x 4 , x 5 ) in the usual linear way. Observe that if we have 
a pentaspherical vector x with nonzero length that corresponds to the origin 
of the underlying cartesian-like coordinate system, a 0(5, C) rotation can 
reduce it to the form x° = (0,0, 0,0, a). If x is of zero length then we can 
reduce it to the form x° = (0, 0, 0, a, —ia) where a is arbitrary Consequently 
if we take the first case and choose a = —i, the pentaspherical coordinates 
(xi, X2, x 3 , x±, —i) map 1-1 to the underlying space via 

X Xi Y %2 Z x 3 

T X4 + 1 ' ^ T £4 + 1' T X4 + I 

where x\ + x% + £§ + x\ = 1, i.e., the unit 3-sphere. For the second case we ob- 
tain, choosing a = | that the pentaspherical coordinates (xi, x 2 , x 3 , 1/2, —i/2) 
map 1-1 to the underlying space via 

x = xi, y = X2, z = X3 

which is tantamount to choosing Euclidean Cartesian coordinates. To con- 
struct superintegrable systems we can invoke what we already know about 
Euclidean space and the three dimensional sphere. There remains the ques- 
tion of equivalence however. 

If we now look at the corresponding problem for the Laplace equation 
there are some modifications necessary. Using the same correspondence for 
pentaspherical coordinates as above, it follows from the work of Bocher [2] 
that if is a solution of the Laplace equation 

(dl + d 2 y + d 2 z +V)1f = 

and if we write ^ = (x 4 + ix 5 ) 1 / 2 $(x, y, z), V = (x 4 + ix 5 ) 2 V(x) the function 
$ satisfies 

(d 2 xi + d 2 X2 + d 2 X3 + d 2 X4 + d 2 X5 + V(x)) $ = 
and has degree of homogeneity —1/2, i.e., 

(xi_d Xl + x 2 d X2 + x 3 d X3 + x 4 d X4 + x 5 d X5 )$ = 

Similarly, the function V has degree of homgeneity —2. Since ^f(x,y,z) = 

*ff( 7^ — , 7^ — , 7^ — ) , the comments made previously about the 

vector x° can be applied. Indeed if we take the case of the three dimensional 
sphere, the Laplace operator acting on the function $ yields the equation 
(A 3 + V — |)$ = 0. Consequently any potential that enables superintegra- 
bility to occur on the complex three sphere also occurs for the corresponding 
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Laplace equation. Similarly, it is clear that if we make the second choice 
for x° we return to superintegrable systems in complex flat space in three 
dimensions. It remains to be shown whether these are all the possibilities. 

To address this, note that the symmetry group acts on the the func- 
tions $ via the normal algebra operators Jy = Xid Xj — Xjd x . and via the 
correct multiplier representation on the functions In the book [2] there 
are mentioned possibilities that there are coordinate systems corresponding 
to S2 x Ei, (the two-sphere x the line). To see how these arise consider 
spheroidal coordinates 

Jpsimp, z = cyj '-{p - l)(p - 1). 



x = Cy/pp cosy?, y 



The corresponding infinitesimal distance is 



ds 2 



dp 2 



p(p 



If we factor out pp from this metric we see that 
ds 2 = c 2 ppds' 2 



c 2 pp 



dp^ 



p 2 (p 



where ds' 2 is clearly a metric corresponding to 5*2 x E\. From what we have 
just done this is conformally equivalent to spheroidal coordinates. However 
it is also conformally equivalent to flat space. Similarly, all such 5*2 x E\ 
systems are conformally equivalent to systems on constant curvature spaces. 



3.4 Theory for 3D classical conformally 2nd order su- 
perintegrable Laplace systems 

Given a classical conformally superintegrable system on a conformally flat 
space we can always find a Cartesian like coordinate system with coordi- 
nates (x,y,z) = (xi,x 2 ,x 3 ) such that the Hamilton- Jacobi (Laplace) equa- 
tion takes the form 

Pl \f + Vl + V-(x) = 0. (26) 
A(x) 

However, this equation is equivalent to the flat space equation 

H=p 2 1 +p 2 2 + P 2 3 + V( X ) = 0, 7(x) = A(x)y(x). (27) 

In particular, the conformal symmetries of (|2"B"1) are identical with the con- 
formal symmetries of ( 1271) . Thus without loss of generality in the classical 
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case, we can assume the manifold is flat space with A = f . (In the quantum 
case a similar result is true but a gauge transformation is required and the 
modification of the potential is dependent on curvature.) 

The operation of inversion in a sphere plays a special role in the classical 
theory. If w(x) is a solution of the Hamilton- Jacobi equation 

Pi + Pi + Pi + ^( x ) = °> Pi = d i u , 

then 

u(x) = Ju(x) = u(^), r 2 = x 2 + y 2 + z 2 (28) 

satisfies 

Pi + Pi + Pi + ^ v (^) = °> Pi = 9iu. 
A second order conformal symmetry 

3 

S = a kj (x)p k Pj + W(x) = C + W, a jk = a kj (29) 
k,j=i 

must satisfy 

{S, H} = 6(x, p)W, b = 6 1 (x)p 1 + 6 2 (x)p 2 + & 3 (x)p 3 , (30) 

for some functions b\, 62, &3- Equating coefficients of monomials in the p's we 
see that the conditions are 

af = 2tf+ai j = 2at + a f < k = ~b i , 
ai + a kl + 0^ = 0, pairwise distinct (31) 

and 

3 

W k = ^a!*V a + afV, k—1,2, 3. (32) 

s=l 

(Here a subscript j on a £m , V or W denotes differentiation with respect to 
Xj.) The requirement that d Xl Wj = d Xj Wi, £ 7^ j leads from f[3"2"j) to the 
second order (conformal) Bertrand-Darboux partial differential equations for 
the potential. 

3 

[V S] a sl - V sl a sj + V 8 {(a s % - (a sj ) e )] + a\% - afV e + (a% - a%)V = 0. 

s=l 

(33) 
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Equations flBTj) are exactly the equations for a second order conformal 
Killing tensor a*- 7 . Thus, necessary and sufficient conditions that S = £ + W 
be a conformal symmetry are that £ is a conformal Killing tensor, W is a 
solution of equations ( 1321) and V satisfies the conformal Bertrand-Darboux 
equations. 

The conformal Killing tensors for flat space are well known, e.g., [51 [T5]. 
The space of conformal Killing tensors is infinite dimensional. It is spanned 
by products of the conformal Killing vectors 

Pi, P'2, P-A, X 3 p 2 ~X 2 p 3 , XiPz-XzPl, X 2 Pl-X 1 p 2 , XlP! + X 2 p 2 + X 3 p 3 , 

(x\ —x\ — x\)pi + 2x\x 3 p 3 + 2x!X 2 p 2 , (xj -x\- x\)p 2 + 2x 2 x 3 p 3 + 2x 2 XiPi, 

(x\ -x\- x 2 2 )p 3 + 2x 3 x 1 p 1 + 2x 3 x 2 p 2 , 

and terms g(x, p)(p 2 +p 2 +P 3 ) where g is an arbitrary function. For a given 
conformal superintegrable system only a finite dimensional space of confor- 
mal tensors occurs. This is for two reasons. First the conformal Bertrand- 
Darboux equations restrict the allowed Killing tensors. Second, on the hy- 
persurface Ti = in phase space all symmetries g(x)7i vanish, so any two 
symmetries differing by g(x)7Y can be identified. 

It is sometimes useful to pass to new variables a 11 , a 24 , a 34 , a 12 , a 13 , a 23 
for the conformal Killing tensor, where a 24 = a 22 — a 11 , a 34 = a 33 — a 11 . Then 
we see that a 24 , a 34 , a 12 , a 13 , a 23 must be polynomials of order < 4. (Thus 
by adding —a u 7i to the second order symmetry we can achieve a 11 = for 
the new tensor with the same action on the hypersurface TC — 0, without 
changing the 5 other variables.) 

For second order conformal superintegrabilty in 3D there must be five 
functionally independent conformal constants of the motion (including the 
Hamiltonian itself). Thus the Hamilton- Jacobi equation admits four addi- 
tional constants of the motion: 

3 

S h=J2 a wPkPi + W {h) = C h + W {h) , h = 1, • • • , 4. (34) 
j,k=i 

We assume that the four functions Sh together with TC are functionally 
independent in the six-dimensional phase space, i.e., that the differentials 
dSh, dTL are linearly independent. (Here the possible V will always be as- 
sumed to form a vector space and we require functional independence for 
each such V and the associated W^ h \ This means that we also require that 
the five quadratic forms Ch, are functionally independent.) We say that 
the functions are weakly functionally independent if dSh, dTi are linearly in- 
dependent for nonzero potentials, but not necessarily for the zero potential. 
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Indeed for now, we will also require that the generating basis is functionally 
linearly independent. 

We can write the system of conformal Bertrand Darboux equations in the 
matrix form Cv = v^Vi + v^V 2 + v^V 3 + v^V, or 






a 12 


a 11 - a 22 


a 31 -a 32 


a 13 





-a 23 


a 21 a 11 — a 


a 32 


-a 32 


-a 13 


a 22 - a 33 a 12 



38 



/ V33 - \ 

v 32 
v 31 



,31 



+ 



a^ 1 + a 
a^ 1 + a 

a 2 — °3 



22 
2 

3.°, 



,32 



a 33 — a 



,38 



,31 



a 22 - a 21 



,ii 



a? 2 — a 



12 



13 
3 
,23 



,22 



,32 



— a 2 / 2 + a 



V 3 + 



33 



,22 
'-21 
,33 
'-31 
,33 
''32 



7 11 

'13 

22 
''23 



v 2 



V. 



(35) 



Corollary 1 Suppose the set {H,Si, ■ ■ ■ ,S^} is functionally linearly inde- 
pendent. Then for general x the 4x5 matrix 



A 



( fl 33 _ 11 

/ "(1) "(1)' 

33 11 

"(2) "(2)' 

33 11 

"(3) U (3)' 

V a 33 - a) 1 



l (4)' 



n 22 - a 11 

a (i) a (i)' 

n 22 - a 11 

"(2) "(2)> 

a 22 - a 11 

"(3) "(3)' 

a 22 - a 11 

"(4) U (4)' 



a 12 
a 12 

"(2)' 

a 12 
"(3)' 

"(4)' 



31 
"(1)' 

31 
"(2)' 

31 
"(3)' 

31 
"(4)' 



a* \ 



32 

(1) 
32 
(2) 
32 

(3) 



n 32 I 

0(4) / 



/ias ranA; ^ ; where the functions aV h Jx) are given by (34\ ). 



There are four sets of equations f[35|) . one for each of the functionally 
independent symmetries (in addition to the Hamiltonian) . We can write 
them as a single matrix equation Bv = b where B is 12x5, b is 12x1 and 





( v 33 - V n 


\ 








V = 


v 12 












\ v 31 


J 



Lemma 1 If the set {H, Si, 

matrix B has rank 5. 



,£4} is functionally linearly independent, the 
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The proof is the same as for the corresponding Helmholtz result in [9J, with 
corrections in the introduction to [llj. 

By choosing a rank 5 minor of B we can solve for v and obtain a solution 
of the form 

V u + A 22 Vi + B 22 V 2 + C 22 V 3 + D 22 V, 
V n + A^V X + B^V 2 + C 3Z V 3 + D 33 V, 

A 12 V X + B 12 V 2 + C 12 V 3 + D l2 V (36) 
A 13 Vi + B 13 V 2 + C 13 V 3 + D 13 V 
A 23 Vi + B 23 V 2 + C 23 V 3 + D 23 V 

If the augmented matrix (B, b) has rank r' > r then there will be r' — r 
additional conditions involving only derivatives less than second order. Here 
the A % \ B lJ , G %3 , D 1 ^ are functions of x that can be calculated explicitly. For 
convenience we take A ij = A j \ B ij = B j \ C ij = C j \ D ij = D ji . 

Suppose now that the superintegrable system is such that r' = r so that 
relations ( 1361) are equivalent to Bv = b. Further, suppose the integrability 
conditions for system (|3"B"1) are satisfied identically. In this case we say that 
the potential is nondegenerate. Otherwise the potential is degenerate. If V is 
nondegenerate then at any point xo, where the B lJ , C l] , D 13 are defined 
and analytic, there is a unique solution V(x) with arbitrarily prescribed val- 
ues of V(xo), Vi(xo), V2(xo), Vs(xo), Vii(xo). The points x are called regular. 
The points of singularity for the A %3 ' , B 13 , C % \ D l i form a manifold of dimen- 
sion < 3. Degenerate potentials depend on fewer parameters. (For example, 
we could have r' = r but the integrability conditions are not satisfied identi- 
cally. Or a first order conformal symmetry might exist and this would imply 
a linear condition on the first derivatives of V alone.) 

Note that for a nondegenerate potential the solution space of fl36l) is ex- 
actly 5-dimensional, i.e. the potential depends on 5 parameters. Degenerate 
potentials depend on < 5 parameters. 



V 22 = 

v 33 = 

v 12 = 

v 13 = 

v 23 = 



3.5 The conformal Stackel transform 

We quickly review the concept of the Stackel transform [S] and extend it 
to conformally superintegrable systems. Suppose we have a second order 
conformal superintegrable system 

H = P \+ P * + f 3 +V(x,y,z) = 0, H = H + V (37) 

and suppose U(x,y,z) is a particular solution of equations flHBT) . nonzero in 
an open set. 
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Theorem 1 The transformed (Helmholtz) system 

H = E, H=(pl + pl + pl)/\ + V (38) 

with potential V(x, y, z) is truly superintegrable, where 

V 

A = XU, V = —, 

PROOF: Let S = ^2 a^PiPj + W = S + W be a second order conformal 
symmetry of TC and Su = Yl a ^PiPj + Wu = <$o + Wu be the special case 
that is in conformal involution with (jp\ + p| +p|)/A + U. Then 

{S, H} = PSo H, {S v , H + U} = p So (H + U), 

and S = S — ^-TC is a corresponding true symmetry of TC. Indeed, 

{S,H} = {S,^}-{Wu^,^} 

This transformation of second order symmetries preserves linear and func- 
tional independence. Thus the transformed system is Helmholtz superinte- 
grable. Q.E.D. 

There is a similar result for first order conformal symmetries £ = ^2 Oj (x)pj 



Corollary 2 Let £ be a first order conformal symmetry of the superinte- 



grable system (Sty and suppose U(x) is a particular solution of equations 
( [HP , nonzero in an open set. Then £ is a true symmetry of the Helmholtz 
superintegrable system [3^1 : {£,71} = 0. 

PROOF: By assumption, {£, TC} = pc(x)TC = PcO^-q + V where pc doesn't 
depend on the momenta pj. Thus, {£, TCq} = pc'Ho, {£, V} = PcV, so a ls° 
{£, U} = p c U Then 

{£, TC} = {£, ^} = H} - U} = Pc (^ - 2<-U) = 0. 

Q.E.D. 

These results show that any second order conformal Laplace superinte- 
grable system admitting a nonconstant potential U can be Stackel trans- 
formed to a Helmholtz superintegrable system. This operation is invertible, 
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although the inverse mapping is not a Stackel transform (it takes true sym- 
metries to conformal symmetries). By choosing all possible special potentials 
U associated with the fixed Laplace system ( 1371) we generate the equivalence 
class of all Helmholtz superintegrable systems ( 1381) obtainable through this 
process. As is easy to check, any two Helmholtz superintegrable systems 
lie in the same equivalence class if and only if they are Stackel equivalent 
in the standard sense. All Helmholtz superintegrable systems are related to 
conformal Laplace systems in this way, so the study of all Helmholtz super- 
integrability on conformally flat manifolds can be reduced to the study of all 
conformal Laplace superintegrable systems on flat space. 

Clearly, the analous results are true in all dimensions n > 2. For n = 2 
they also extend to operator Laplace equtions without change; for n = 3 they 
also extend but, as discussed above, the potentials must be modified under 
the Stackel transform to take scalar curvature into account. 



Example 4 Consider the degenerate Laplace system (d)] in three variables 
x,y,z. If we choose U = 1/z 2 and perform the corresponding Stackel trans- 
form we obtain a Helmholtz superintegrable system on the 3-sphere with two- 
parameter potential. This follows immediately from the fact that the metric 
(dx 2 + dy 2 + dz 2 ) / z 2 corresponds to a space with nonzero constant cuvature. 



3.6 The integrability conditions for the potential 



To determine the integrability conditions we first introduce the dependent 
variables W {0 \ W w = V u W (2) = V 2 , = V 3 , W w = V n , the vector 



w 



.tr 



(W^,W^\W^,W^,W ( ^), 



(39) 



and the matrices 



A (D 



A (2) 
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A 12 
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D 13 


A 13 


£13 
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B u 
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12 _ ^22 
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12 B 


12 (J 


12 
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22 B 
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23 B 
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24 A 


24 B 


24 (J 


24 


A 12 J 
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A^ = 



/ o 

D 13 
D 23 

\ D 34 





A 13 B 13 

A 23 B 23 

A 33 B 33 

A 34 B 34 



1 

C 13 
C 23 

c 33 



o \ 




1 



(42) 



C 34 A 13 J 



where 

A 14 -- 

B 14 -- 

C 14 = 

D 14 = 

A 24 = 

C 24 -- 

A 34 -- 

C 34 -- 



A\ 2 
B] 2 
C\ 2 
D\ 2 



A 22 + B 12 A 22 + A 12 A 12 - B 22 A 12 - C 22 A 13 + C 12 A 23 - D 22 , 
Bf + A 12 B 12 - C 22 B 13 + C 12 B 23 + D 12 , 
C 22 + B 12 C 22 + A 12 C 12 - B 22 C 12 - C 22 C 13 + C 12 C 23 
Df + A 12 D 12 - B 22 D 12 + B 12 D 22 + C 12 D 23 - C 22 D 13 , 
A 12 + B 12 A 12 + C 12 A 13 + D 12 , B 24 = Bf + B 12 B 12 + C 12 B 13 , 

12 , D 12^12 , C 12 C 13 j jD 24 = jB 12 jD 12 + c ,1 2jD 13 + jD 12 ^ 



Cf + B V2 C V 
A\ 3 + B 13 A 12 + C 13 A 13 + D 13 , 

C 13 + ^13^12 + C 13 C <13 ) D 34 



B 



34 



13 D 13 



5^ + + C Li B 



D 13 + B 13 D 12 + C 13 D 13 . 
Then the integrability conditions for the system 

d x .w = A®w j = 1,2,3, 

must hold. They are 

Af - Af = A (i) A (i) - A (j) A w = [A ( f A (i) ] . 

For convenience in the arguments to follow we set 

U 1 = Af - Af - [A( 2 \ A< 3 >], U 2 = Af - Af - [A®, A«], 

W 3 = Af -Af -[A«,A( 2 )], 
so that the identities are 

U 1 = U 2 = U 3 = 0. 
The simplest of these, i.e., those that don't involve derivatives, are 

C 13 



(44) 



(45) 



(46) 



(47) 



A 23 
B 23 - A 13 



B 13 

C 22 , 



C 12 , B 12 
A 12 + B 33 . 



A 22 + A 33 
C 23 . 



Thus we can write all of the A % \ B % \ C u in terms of the 10 functions 
A 12 , A 13 , A 22 , A 23 , A 33 , B 12 , B 22 , B 33 , C 33 . 



(48) 



(49) 
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Also, identities (H7|) enable us to express each of D 33 , D 23 , D 22 , D 13 , D 12 as a 
polynomial in the A*- 7 , B % \ C u and their first derivatives for i, j < 3. 

As an example, the nondegenerate potential IHM satisfies the canonical 
equations ( |36i) with 



^33 = 3 + 12x(z 2 -x 2 ) ^ ^3 = I2y(z 2 -x 2 ) ^ ^ = _3 + 12^ 2 -x 2 ) ^ 
x 1 — r 4 1 — r 4 2; 1 — r 4 

(50) 

r33 24(z 2 -x 2 ) ^ 22 3 | l&fo'-s 2 ) ^ 22 3 | 12y(y 2 -a; 2 ) ^ 
1 — r 4 ' a; 1 — r 4 y 1 — r 4 

g22 _ 12z(y 2 - a; 2 ) d22 _2% 2 -x 2 ) 



A 23 
A 13 
A 12 



1 - r 4 1-r 4 

1 2xyz 
1-r 4 ' 
12x 2 ^ 



1 -r 4 ' 
12x 2 y 



s 23 = 


Yly 2 z 


C 23 = 


Ylyz 2 


D 23 = 


24yz 


1 -r 4 ' 


1-r 4 ' 


1-r 4 


B 13 = 


12:q/z 


C 13 = 


12xz 2 


D 13 = 


24xz 




1-r 4 ' 


1-r 4 ' 




1-r 4 


B 12 = 


Ylxy 2 


C 12 = 


Ylxyz 


D 12 = 


24xy 


1-r 4 ' 


1-r 4 ' 


1-r 4 



l_ r 4> 

where r 2 = x 2 + y 2 + z 2 . The nondegenerate potential (1221) satisfies a simialr 
set of canonical equations but the expressions for the terms analogous to (I5U|) 
are somwhat more complicated and would take two pages to list. 



3.7 Integrability conditions for the symmetries 

Since (as we assume) the potential is nondegenerate, at any regular point 
Xq, V and the first derivatives V±, V2, V3 can be chosen arbitrarily. Thus the 
coefficients of V and Vj on both sides of equation (J3"5"|) must be equal. From 
this, we obtain the relations 

2a 13 = a 13 D 33 + (a 11 - a 33 )D 13 + a 12 D 23 - a 23 D 12 , 

_3a 31 = -a 12 A 23 + (a 33 - a n )A 13 + a 23 A 12 - a 13 A 33 , 

a 12 _ a 32 = _ a l 2j g23 + (fl 33 _ a ll )jB 13 + fl 2 3 £l2 _ 

3a 13 = -a l2 C 23 + (a 33 - a n )C 13 + a 23 C 12 - a 13 C 33 , 

with 8 analogous relations from the other two Bertrand-Darboux equations. 
Using these 12 relations and equations (I3T1) we can solve for all of the first 
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partial derivatives a{ for j ^ k and a" — a 3 / to obtain 

3a} 2 = a 12 A 22 -(a 22 -a n )A 12 -a 23 A 13 + a 13 A 23 , (51) 

3(a n -a 22 ) 2 = 2[-a 12 A 22 + (a 22 -a u )A 12 + a 23 A 13 -a 13 A 23 }, 

3af = - a 12 C 23 + (a m -a n )C 13 + a 23 C 12 -a 13 C 33 

3(a 33 -a 11 ) 1 = 2[a 12 C 23 -(a 33 -a n )C 13 -a 23 C 12 + a 13 C 33 ], 

3a 23 = a 23 (B 33 -B 22 )-(a 33 -a 22 )B 23 -a 13 B 12 + a 12 B 13 , 

3(a 22 -a n ) 3 = 2[-a 23 (A 12 + J B 33 - J B 22 ) + (a 33 -a 22 ) J B 23 + a 13 (S 12 + A 33 ) 
+a 12 (A 23 - B 13 ) + (a 11 - a 33 )A 13 ] , 

3a} 3 = -a 23 A 12 + (a n -a 33 )A 13 + a 13 A 33 + a 12 A 23 , 

3(a 33 -a n ) 3 = 2 [-a 23 A 12 + (a 11 - a 33 )A 13 + a 13 A 33 + a 12 A 23 } , 

3(a 33 -a n ) 2 = 2[a 13 (A 23 -C 12 ) + (a 22 -a 33 )C 23 + (a 11 -a 22 )A 12 
+a 12 (A 22 + C 13 ) - a 23 (A 13 + C 22 - C 33 )] , 

3a 23 = a 13 C 12 -(a 22 -a 33 )C 23 -a 12 C 13 -a 23 (C 33 -C 22 ), 

3al 2 = -a 13 B 23 + (a 22 -a n )B 12 -a 12 B 22 + a 23 B 13 , 

3(a 22 -a 11 ) 1 = 2[a 13 S 23 -(a 22 -a 11 ) J B 12 + a 12 J B 22 -a 23 J B 13 ], 

3af = a 12 (B 23 + C 22 ) + a u (B 13 + C 12 )-a 22 C 12 -a 33 B 13 

+ a 13 (B 33 + C 23 )-a 23 (C 13 + B 12 ), 

3af = a 12 (-2B 23 + C 22 ) + a u (C 12 -2B 13 )-a 22 C 12 + 2a 33 B 13 

+ a 13 (-2B 33 + C 23 ) + a 23 (-C 13 + 2B 12 ), 

3al 3 = a 12 (B 23 -2C 22 ) + a 1 \B 13 -2C 12 ) + 2a 22 C 12 -a 33 B 13 

+ a 13 (B 33 -2C 23 ) + a 23 (2C 13 -B 12 ). 

There are several conditions left over. These are the obstructions 

a 12 (C ,22 _ £23 + A 13J + (a ll _ a 22 )(C 12 _ ^23) + _ a 33 )(A 23 _ ^13) (53) 

+a 13 (C 23 - B 23 - A 12 ) + a 23 (B 12 - C 13 + A 33 - A 22 ) = 

and 



2a} 2 


= a 12 L> 22 + l 


[a 11 - 


- a 22 )D 12 


+ a 13 D 23 -a 23 D 13 , 


(53) 


2a} 3 


= a 13 L> 33 + l 


y 1 - 


• a 33 )D 13 


+ a 12 D 23 -a 23 D 12 , 


(54) 


2a 23 


= a 23 (D 33 - 


D 22 ) 


+ (a 22 - 


a 33 )D 23 + a 12 D 13 -a 13 D 12 . 


(55) 



It follows directly from conditions ( HBl) for a nondegenerate potential that 
obstruction (1521) is satisfied identically. 
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3.8 Nondegenerate potentials with 6 linearly indepen- 
dent second order conformal symmetries 

Suppose we have a superintegrable system with nondegenerate potential and 
5 functionally independent second order symmetries Si,-- - ,Si,H. We al- 
ready know that any second order superintegrable Helmholtz system with 
nondegenerate potential on any conformally flat space will lead to such a 
Laplace system, in fact a system with 6 linearly independent second or- 
der symmetries. Now we will demonstrate that, conversely, every Laplace 
system with nondegenerate potential leads to a Stackel equivalence class of 
Helmholtz superintegrable systems on conformally flat manifolds. Let 

n = P 2 1 + P 2 2 + P l + v( a , x ) = o 

be a nondegenerate Laplace system, where the parameters are denoted by 
the vector a = (ai, • • • ,05). Let A(x) = V^(ao,x) be a special case of this 
potential with fixed parameters. With a suitable linear transformation in 
parameter space we can always assume 

V(a, x) = A(x)u(e, x) — a 5 A(x), a=(e,a 5 ). (56) 

By assumption, both V and its special case A satisfy the canonical equations 
for the potential. Substituting V = XU in these equations we see that 
the Stackel transformed potential U satisfies the canonical equations 

U22 = U u + A 22 Ui + B 22 U 2 + C 22 U 3 , 

C/33 = U u + A 33 Ui + B 33 U 2 + C 33 U 3 , 

U 12 = A l2 Ui + B l2 U 2 + C 12 U 3 , (57) 

Ui 3 = A 13 U! + B 13 U 2 + C 13 U 3 , 

U 23 = A^ + B^U. + C 23 ^, 

characteristic of a Helmholtz system with nondegenerate potential, where 

A 33 = A 33 + 2^, B 33 = B 33 , C 33 = C 33 - 2^, A 22 = A 22 + 2^, 
A A A 

B 22 = B 22 — 2— C 22 = C 22 A 12 = A 12 — — B 12 = B 12 — — C 12 = C 12 



A 2 
2— 

A ' 


c 22 = 


fj22 j^12 _ p\2 


A 2 
A ' 


i 13 


= A 13 


_ b 13 = B 13 
A 


, c 


A 23 


= A 23 


B 23 — b 23 — — 


, c 



A 



j: i / < : 3 _ ^± 

A ' 

23 _ /^f23 ^2 



28 



In particular, the coefficients of U vanish in all these equations. The corre- 
sponding Stackel transformed system is 

H^ = pl+P } +P h u = Pl+P } +P h v-a 5 = 0. (58) 
A A A 

This appears to be a nondegenerate Helmholtz superintegrable system with 
energy 05. However, we know only that <Si, • • • , 1S5 are conformal symmetries 
of this system, not the required true symmetries. We need to exhibit true 
symmetries. Recalling the conclusion of Theorem [TJ if S is a second order 
conformal symmetry of TC = then S = S — is a true symmetry of 

(158j) . Thus the new conformal symmetries 

W x 

<Sh = <Sh h = 1, - ■ ■ ,5, 

are actually true symmetries, i.e., in involution with the Hamiltonian (p\ + 
P2 + P3) A + t '( e ; x )- Note that the symmetries Sh, Sh agree on the hypersur- 
face Ti — and {Sh,S t } vanishes on the hypersurface if and only if {Sh, S t } 
vanishes. 

We see that a Laplace superintegrable system with nondegenerate po- 
tential and a guaranteed 5 functionally independent second order conformal 
symmetries is equivalent via a Stackel transform to a Helmholtz system with 
nondegenerate potential and 5 functionally independent second order true 
symmetries. However the 5 ==>- 6 Theorem in [9] shows that such a Helmholtz 
system actually admits 6 linearly independent true symmetries. This extra 
symmetry must correspond to a conformal symmetry of the original Laplace 
system. Thus the Laplace system must admit 5 conformal symmetries, in 
addition to the Hamiltonian. 

Theorem 2 (4 ==>• 5) A Laplace second order superintegrable system with 
functionally independent generators Si,-- - ,5 4 ,7i, and nondegenerate po- 
tential admits a 5th conformal second order symmetry S§ such that the set 
{Si, ■ ■ ■ ,Ss} is linearly independent on the hypersurface TC — 0. 

Theorem 3 There is a one-to-one relationship between flat space Laplace 
systems with nondegenerate potential and Stackel equivalence classes of su- 
perintegrable Helmholtz systems with nondegenerate potential on conformally 
flat spaces. 

For such a Laplace system the integrability conditions for the potential (j47p 
and the integrability conditions for the symmetries (1511) . some 150 equations 
in all, are satisfied identically. Furthermore each of the obstruction equations 
(1531) . ( 1541) . ( 1551) must also be satisfied identically. 
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By a straightforward but lengthy calculation with MAPLE we established 
the following from these equations: 

1. We found 30 equations expressing the 30 partial derivatives diF^ k as 
quadratic polynomials in the 10 functions A 12 , A 13 , A 22 , A 23 , A 33 , B 12 , 
B 22 , B 23 , B 33 , C 33 . Here F J is any one of these functions. For example, 
three of these equations are 

q A 12 = 1^23^13 + A 23 B 23 ^33^33 _ 1^22 _ 1^12^12 

3 3 3 

— B 33 A 22 — A 23 C 33 

dA 12 = I (A 23 )2 I (A 12 )2 I (5 12 )2 I (A 33 ) 2_I c ,33 A 13_I (jB 33 ) 2 

y 2 6 2 6 6 6 V ; 

_ 1^22^33 _ ^^33^12 , 1^33^12 , 1^33^22 _ ^/^>23\2 , ^^33^23 

6 3 3 6 6 l ; 6 

d A 12 = i A 23 A 33 + 2 1 2 ^23 1^13^12 _ 1^22 

3 3 3 3 

2. We found quadratic expressions for each of the terms D % K Indeed: 

D 12 = 2 ( _ a i 2jB i 2 + a2 3 jB 23 + jB 33 A 33_ jB 33 A 22_ a 23 C 33 )) 

3 

D 13 = 2 ( _ A 13 S 1 2 _ S 2 2A 23_ 5 23 A 33 + 5 23 A 22 + A 23 5 33 + A 12 A 23 ); 

3 

D 22 = 2 ^ 12)2 + ^12^22 _ £,22^12 + ^33^13 + 2 B 33 A 12 - 2A 33 B 12 - C 33 
_ B K B 22 + (£ >2 3) 2 + (j g33 )2 _ (£12)2 _ ^33)2 + ^33^22 _ £23^13) 

D 23 = -(A 23 A 33 + B 12 A 23 - B 23 A 12 - A 13 B 33 ), 
3 

D 33 = -(B 33 A 12 — (B 12 ) 2 — A 33 B 12 — C 33 B 23 — B 33 B 22 + (B 23 ) 2 + (£? 33 ) 2 ) 
3 

3. There are no other consequences of these identities. 

There are no polynomial identities that the 10 functions must obey. What 
is amazing is that the integrability conditions for the expressions diF^ h are 
satisfied identically! 

Theorem 4 Suppose the integrability conditions (ffi) and the integrability 
conditions for the symmetries (f57]). are satisfied identically. Then equa- 
tions / fJPj) hold and the integrability conditions de(diF^ k ) = d^d^F^ ) are 
satisfied identically. Thus at a regular point x we can choose a 10-tuple 
c = (ci, • • ■ , cio) arbitrarily and there will exist one and only one superinte- 
grable system such that 

(A 12 (x ), A 13 (x ), • ■ ■ , 5 33 (x ), C 33 (x )) = c. 
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3.9 The classification problem for nondegenerate po- 
tentials 

Theorem H] provides the basis for classifying all manifolds for which Helmholtz 
superintegrable systems exist, and all nondegenerate potentials on these man- 
ifolds, a program that is not yet complete. Indeed, note that the conformal 
group, with connected component isomorphic to 5*0(5, C ), acts naturally 
on the 10-tuple c. Suppose we have a conformal superintegrable system 

p\+pl+pl + V{^) = 0. (60) 

If g : x — > x' = xg is an element of the conformal group (considered as a 
transformation group) then g acts on functions /(x) via operators T(g) such 
that T(p)/(x) = /(xp). Then T(g 1 g 2 ) = T(g 1 )T(g 2 ) so we get a representa- 
tion. Under this action p\ +p 2 + p\ — > p' 2 +p 2 2 +p'z = c(x, g) (pi 2 +p 2 2 +P3 2 ) 
where c(x, g) is the conformality factor. Thus the conformally superinte- 
grable system transforms to another conformally superintegrable system 

pi 2 + P2 2 + p 3 2 + H*) = 0, Hx) = -r—, v ^ (6i) 

The conformal symmetries transform in an obvious manner. System f loT)]) is 
uniquely characterized by the values of the 10-tuple c at the regular point 
xo. There is an induced action g : c — > c' = eg of the conformal group 
on 10-tuples such that system ( 16T1) is uniquely determined by the values 
c g at the point x g. Since the equations determining this action locally 
on 10-tuples are autonomous for Euclidean actions, we can mostly ignore 
the starting point xo and focus just on the map c — > eg. However, for 
general conformal actions we have to consider the map on the 13-dimensional 
manifold of points (x, c). Thus (xo,c) — > (xog,cg). Clearly, all conformally 
superintegrable systems related by elements of the conformal group have 
essentially the same structure and should be identified. Thus we say that 
two superintegrable systems are equivalent if and only if they are on the 
same orbit under the action (x, c) — > (xg,cg). It appears that we might be 
able to determine a solution in each equivalence class and thus, indirectly, 
find all 3D nondegenerate Helmholtz superintegrable systems, including all 
those on nonflat spaces. 

To see how the classification might proceed, recall from [0] that every 
Helmholtz superintegrable system with nondegenerate potential is Stackel 
equivalent to a system on either flat space, or the complex 3-sphere, or both. 
Thus to classify all possible systems we need only to classify the constant 
curvature space systems. From the point of view of this paper, the flat 
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space systems are just those for which Z7 y = in equations (|59|) . Thus 
1^ = ■ ■ ■ = 1^ = at each regular point, where 

jW = -A 22 B 23 + B 23 A 33 + B 12 A 13 + A 23 B 22 - A 12 A 23 - A 23 B 33 (62) 
j(b) _ (^ 33 ) 2 _|_ B 12 A 33 — A 33 A 22 — A 12 B 33 — A 13 C 33 + A 12 B 22 
-B 12 A 22 + A 13 B 23 -(A 12 ) 2 

j(c) = B 23 C 33 + ^12^33 + (£12)2 + ^22^33 _ (j g3 3) 2 _ ^12^33 _ (£23)2 

j(d) = _ j gl2^23 _ yp^23 + ^13 £33 + ^12 £23 

j(e) _ _£23^23 _j_ ^33^23 _j_ ^22^33 _ ^33^33 _|_ £12^12 

In [12] it was shown via Grobner basis methods that these conditions further 
imply 1^) = where 

/(/) = A 13 C 33 + 2A 13 B 23 + B 22 B 33 -(B 33 ) 2 + A 33 A 22 -{A 33 ) 2 (63) 
+ 2A 12 B 22 + (A 12 ) 2 - 2B 12 A 22 + (B 12 ) 2 + B 23 C 33 - (B 23 ) 2 - 3(A 23 ) 2 , 

and that any system satisfying conditions 1^ = ■ ■ ■ = = at some 
regular point uniquely defines a flat space superintegrable system. From this 
we were able to classify all flat space Helmholtz superintegrable systems. To 
compare our present notation to the results of [9], consider Laplace super- 
integrable systems with nondegenerate potential that are Stackel equivalent 
to a Helmholtz superintegrable system on the complex 3-sphere. It is always 
possible to choose coordinates x, y, z on the 3-sphere such that the metric 
takes the form X(x,y, z) = 4/(1 + x 2 + y 2 + z 2 ) 2 as in (fTB"]) . Thus the non- 
degenerate potential for the associated flat space Laplace system must have 
V = X(x,y,z) as a particular instance. Substituting this requirement in 
equations (1361) we can solve for the functions to get 

D 12 = -h^ = G xy + G x G y -A 12 G x -B 12 G y -C 12 G z , (64) 

-D 13 = ~\l {a) = G xz + G x G z -A 13 G x -B 13 G y -C 13 G z , 

D 22 = -h^ + h^=G yy + G 2 y -G xx -G 2 x -A 22 G x -B 22 G y -C 22 G z , 
o 

r-)23 j(d) /~i \ r 1 r 1 A^C* C?23/^( ri2"i>ri 

D — — -1 — Ur yz -+- (J y (j z — 1\ Kj x — Li Kj y — <_/ Ur z , 

D 33 = -h^ = G zz + G 2 z -G xx -G 2 x -A 33 G x -B 33 G y -C 33 G z , 

where G(x, y, z) = In A. Equations (I5T?]) and (I6"4l) should play an impor- 
tant role in the classification of all Helmholtz superintegrable systems on the 
complex 3-sphere with nondegenerate potential. 
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